Late stages of coarsening in model C 
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We present a comprehensive picture of (non-critical) domain growth in model C systems where 
a non-conserved scalar order parameter is coupled to a conserved concentration field. For quenches 
into the region where the ordered and disordered phases coexist, we confirm earlier partial numerical 
results and find a growth exponent z = 3. For quenches into the ordered region, we confirm the 
theoretical prediction z — 2. Finally we discuss the implications of our results for domain growth in 
the microcanonical (/)^-model and we offer some criticism of the work of Somoza and Sagui on the 
morphology and wetting properties of domains. 



I. INTRODUCTION 



Coarsening has been a object of intensive study over 
the last decades not only because of its experimental rel- 
evance but also because it constitutes what is maybe 
the simplest case of out-of-equilibrium dynamics. For 
a scalar order parameter, one generally distinguishes the 
non-conserved case, e.g. the Ising model with Glauber 
dynamics (model A according to the widely-used classi- 
fication of [1]), from the conserved case, e.g. the Ising 
model with Kawasaki dynamics (model B). It is now 
well-established that in these cases phase ordering is 
characterized by a single lengthscale L growing alge- 
braically with time {L ~ t^^^) with the growth expo- 
nent z taking two different values, z = 2 for the non- 
conserved/model A case, and z = 3 for conserved order 
parameter systems like model B [2]. Here we are in- 
terested in the more complicated and therefore less well 
understood case of a non-conserved order parameter cou- 
pled to a conserved concentration (so-called model C). 
Examples of this situation can be found in various phys- 
ical systems, e.g., intermetallic alloys (see [3] and refer- 
ences therein), adsorbed layers on solid substrates [4] and 
supercooled liquids [5]. In these systems, the decomposi- 
tion process (described by the conserved concentration c) 
and the ordering process (described by the non-conserved 
order parameter i/j) are coupled. Let us consider an A — 5 
alloy on a square lattice in order to illustrate this. A 
state in which all A atoms are surrounded by B atoms 
is energetically favorable. The ordered state thus con- 
sists of two symmetric sublattices, and we can define an 
order parameter ip as half of the difference between the A- 
concentration in each sublattice. In this way, = 1 when 
all the A atoms are on the one sublattice and ip = —1 
when they are on the other. At high temperature, a dis- 
ordered state tjj = arises. It is now easy to realize that 
for asymmetric initial conditions (i.e. an unequal amount 
of A and B atoms) the system will not be able to com- 
pletely order (strictly speaking, this is only true at low- 
enough temperature). Hence, as opposed to model A, 
the disordered phase can coexist with the two ordered 
phases. On a typical equilibrium phase diagram in the 



concentration-temperature (c-T) plane (Fig. 1), one can 
thus distinguish, apart from a disordered region and an 
ordered region, a coexistence region. The dashed line 
separating the ordered and disordered regions marks a 
second-order phase transition. In the spinodal region in- 
side the coexistence region (dotted line), the three phases 
are thermodynamically unstable. 
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FIG. 1. Typical equilibrium phase-diagram of a model G in 
the (c, T) plane. The solid line is the coexistence curve, the 
dashed curve is the order-disorder transition line. The dotted 
line within the coexistence region is the spinodal curve inside 
which the three phases are unstable. Four different types of 
quenches (Ql, Q2, Q3, and Q4) are shown. 

Models have been proposed to account for various as- 
pects of the morphology and of the kinetics of the ex- 
perimental systems (see for instance [3] and references 
therein). From the more theoretical point of view of uni- 
versality issues, the situation is not quite satisfactory. 
For instance, the critical exponents, and in particular 
the dynamic critical exponent, are still debated [6-8]. A 
renormalization group analysis turns out to be more del- 
icate than in the case of model A [9,10]. Our goal here is 
to clarify the a priori simpler problem of domain growth 
below criticality, when the system is quenched down from 
a high-temperature state. Notable but partial results, 
somewhat scattered in the literature, have been obtained 
in the past. For quenches into the spinodal region with 
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droplet morphology (quench Q2 of Fig. 1) San Miguel et 
al. [11] have predicted the model B exponent z = 3. Nu- 
merical simulations in the context of a Oono-Puri "cell 
model" have been found to be consistent with this pre- 
diction [12,13]. On the other hand, Elder et al. [14] have 
predicted z = 2 for quenches above the tricritical temper- 
ature, i.e. in the ordered region (quench Q4). To the best 
of our knowledge, this has not been verified numerically. 

Our goal here is to give a complete picture of (non- 
critical) domain growth in model C, considering, within 
a single system introduced in Section II, all four possible 
types of quenches illustrated in Fig. 1. This is done in 
Section III. In Section IV, in the sake of comprehensive- 
ness, we come back to the two following unsettled issues 
discussed recently in works about model C systems. The 
microcanonical (j)^ model [15,16,7], is a type of model C 
since the order parameter is coupled to the (conserved) 
energy. Zheng has suggested in a recent paper [17] that 
domain growth is characterized by a non-trivial value of z 
(2 < z :^ 2.65 < 3). A more careful study by us showed 
that the data are in fact consistent with the model A 
exponent z = 2 [18]. Here we detail to which phase 
of model C the microcanonical (j)^ model belongs. The 
morphology of domains and the related "wetting" issues 
have also been a point of contention in the past. In ex- 
periments, it has been observed that neighboring ordered 
domains do not merge [19]. A possible explanation pro- 
posed in [20] is that the domains are different variants of 
the same ordered structure. The simulations of [3] seem 
to indicate that ordered domains do not join but "stay 
separated by narrow channels of the disordered phase": 
the antiphase boundaries appear to be wetted by the dis- 
order phase. But Somoza and Sagui [22] have found on 
the contrary that inside the coexistence region the two 
ordered phases may be in direct contact. We revisit their 
work and resolve the controversy. A summary of our re- 
sults is given in Section V. 



II. THE MODEL 

We choose one of the simplest versions of model C 
which can be written as follows: 



~dt 



dt ' (5c 



(1) 

(2) 



Here and Tc are kinetic coefficients, rj and rj' rep- 
resent thermal noise and F [1/^,0] is the Ginzburg-Landau 
free energy functional which takes generally the following 
form: 
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(Vc)2 



(3) 



where Kc and are diffusion constants. The function 
f{c,tlj) should satisfy a few constraints. Firstly it has to 



be symmetric in ip: f{c^ip) = f{c,—ip). It should also 
allow for the coexistence of a disordered phase (cn,'0n) 
with the ordered phases (cg, ±V^s) [21]. These points cor- 
respond to minima in the free energy landscape. A pos- 
sible choice is [22,23] 



/(c,tA) = r(i-r)' + «(r-c) 



(4) 



where a governs the coupling between the two fields. The 
minima are then (cn,'0n) = (0^0) corresponding to the 
disordered phase and (cg, ?/^s) = (1, ±1) corresponding to 
the ordered phase. 

It is instructive to analyze the model in terms of its uni- 
form fixed points (ccV^o)- They lie on the curves where 
fip=0 which are the c-axis i/jq = and the polynomial 
^0 = ~ ^^o)(l - V^o) + ^0- Their stability is deter- 
mined by the signs of /^^ and {fccUi; - f%) [11,23]. 

Fixed points on the c-axis change stability at the so- 
called ordering spinodal c\\ they are stable for c < 
ci = The point on the second curve where stability 
changes is the conditional spinodal c^r- Stability occurs 
when c > Ccr 
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FIG. 2. (c,'0)-plane. Solid lines: stable uniform fixed 
points. Dotted lines: unstable uniform fixed points. The 
minima of / are indicated by filled circles. Grey lines de- 
limit the coexistence domain. Dashed grey lines delimit the 
spinoidal region 

Consider now the evolution of a two-dimensional sys- 
tem starting from random initial conditions disordered 
around = with some fixed concentration c^. We 
briefiy present the expected follow-up of quenches at var- 
ious Ci values before describing our corresponding numer- 
ical results in the next section. 

— For Ci < Cn (disordered phase), the system evolves 
towards the stable point (ci,0). 

— For Cn < Ci < Cs, the three phases (0,0) and (1,±1) 
can in principle coexist. 

— For Cn < Ci < c\ (quench Ql), however, (ci,0) is still 
stable and it will be reached if initial fiuctuations are 
small enough. 
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— Inside the region ci < q < Ca- (quench Q2), the sys- 
tem necessarily undergoes spinodal decomposition. For 
Ca < Ci < Cs (quench Q3), even though the ordered fixed 
points are now stable, nothing is changed since the initial 
conditions (chosen around = 0) lead to the formation 
of disordered domains. 

— For initial concentrations q > c (quench Q4) , the sys- 
tem will not phase separate but will order, so we identify 
this region with the ordering regime. 



III. NUMERICAL RESULTS 

In order to investigate the late stages of domain 
growth, Eqs. (1-2) are now studied numerically. We 
perform quenches from the high-temperature disordered 
phase to T = 0. (The noise term in Eqs. (1-2) can be set 
equal to zero, since the growth exponent is expected to 
be determined by the "T = fixed point" [2].) In prac- 
tice, initial conditions are randomly distributed around 
(ci,0). More specifically, we choose ip G [—0.1 : 0.1] and 
cG [q - 0.05 :ci + 0.05]. 

Equations (1-2) are numerically solved on a two- 
dimensional grid. For the time integration we use Euler's 
method and we approximate the Laplacian by: 
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Ax2 



(5) 



where V/ and Vf are the sets of nearest and next-nearest 
neighbors of site i. 

In most regimes studied below. Ax = 1.0 and At = 
0.01 allow for a smooth representation of both fields. 
Without loss of generality, we set: Kc = = 1 and a = 
4. We record the typical domain size L of both the order 
parameter and the concentration fields, as determined by 
the mid-height value of Cc,^(r, t), the normalized two- 
point correlation function calculated for simplicity along 
the principal axes of the lattice using the the reduced 
"spin" variables = sign('0) and dc = sign(c — q): 



shows that the spinodal lines do not have a strict mean- 
ing in the presence of fiuct nations, a well-known fact (see 
for instance [24]). 



When decomposition occurs, the ordered phase is the 
minority phase, and, similarly to what is observed in 
model B, droplets of either ordered phases form. A sim- 
ilar configuration is shown in Fig. 4. The growth of or- 
dered droplets is constrained by the conservation law. 
The smallest droplet disintegrates to the benefit of the 
larger ones. We thus expect a Lifshitz-Slyozov growth 
law L ~ t^/^. In our simulations, this behavior is not 
easily observed, due to the presence of a long transient 
during which the ordered domains "nucleate" out of the 
met ast able state. In particular, the growth of the typical 
lengthscales Lc,^ does not show the expected scaling al- 
though the trend at large times is good. In this regime, 
though, the droplet morphology of the ordred domains 
can be exploited to provide a more direct measurement 
of growth. In Fig. 3, we show the evolution of the total 
number of the ordered droplet, s n, and of the average 
droplet volume V. These quantities scale respectively as 
and t^/^. Due to a long crossover, the expected 
behavior is only clearly seen when plotting V and 1/n 
against t^/^ (Fig.Sb). 
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100 200 300 400 
3. Domain growth at a = 0.1. (a) V and l/n (dot- 
ted lines) vs. t in log- log plot. The dashed lines indicate the 
expected asymptotic behavior V,l/n t^^^ (b) V and l/n vs. 
t^/^ (1/ n is multiplied with some factor so that it can diplayed 
together with V). System size is 2048 x 2048. 



A. Quenches into the coexistence regime: droplet 
regime 

Quenches within the disordered regime being uninter- 
esting, we first discuss Ql quenches, i.e. quenches into 
the coexistence region, but left of the spinodal (c^ < 
c < ci). Here the system possesses a stable fixed point 
{ci , 0) but will be able to locally reach the ordered fixed 
points if initial fiuctuations are large enough. The precise 
dependence on initial conditions is in fact a non-trivial 
problem that we have not studied systematically. This 



The same scenario is observed for Q2 quenches pro- 
vided the ordered phase be the minority phase, i.e. for 
ci < Ci < ^^2^^ =0.5. A typical snapshot is shown in 
Fig. 3. We have also verified that indeed a growth expo- 
nent 2; = 3 is observed. Now the crossover is shorter and 
we see that the local slope of log Lc,^ vs logt approaches 
I at late times (see Fig. Sab). The same conclusion is 
reached studying the total number of domains and their 
average volume (Fig. 5cd). 
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FIG. 4. Snapshots of coarsening in droplet regime 
(a = 0.25) for a system of 128 x 128 at t = 200. High val- 
ues of the field correspond to dark areas, (a) concentration 
field: white corresponds to the disordered phase c = 0, black 
to the ordered phase c = 1, (b) order parameter field: grey 
corresponds to the ordered phases '0 = 0, white and black 
correspond to the ordered phases ^ = ±1. 
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FIG. 5. Domain growth in the droplet regime (ci — 0.25) 
(a) Lc vs i (b) vs t (c) V and 1/n (dotted lines) vs t (d) 
V and 1/n vs t^l'^ . (a),(b) and (c) are on logarithmic scale, 
(d) is in linear scale. The dashed lines indicate the expected 
scaling / exponent values. Insets: local slopes calculated by a 
running average of the derivated signal over a window of fixed 
size in the corresponding variable (i.e. log^g ^ '^^^ here). 
System size is 4096 x 4096. 



B. Quenches into the coexistence regime: 
interconnected regime 

Due to the coupling to the non-conserved order param- 
eter, the symmetry between high and low concentration 
phases, which exists in model B, is broken. For Ci > 1/2 



the minority phase (now the low concentration phase) 
does not form droplets but is instead concentrated along 
lines which correspond to interfaces of the non-conserved 
field. In the context of alloys one speaks of "wetting 
of antiphase boundaries" . Looking at the non-conserved 
field, we first observe model A like coarsening which then 
slows down, since the inter facial regions separating or- 
dered domains get wider (made of the disordered phase, 
their volume is conserved). We thus again encounter the 
situation where the ordering is constrained by the con- 
servation law, but with the difference that the disordered 
phase has a particular morphology. It is therefore not ob- 
vious that the growth exponent will be 2; = 3. The deriva- 
tion of [11] is only valid in the case of isolated droplets. 
However, the Lifshitz-Slyozov law is known to apply more 
generally, provided that 1) coarsening is driven by surface 
tension, 2) transport is by diffusion through the bulk and 
3) the length scale that describes the coarsening process 
is the only relevant length scale in the system [25]. It is 
easy to convince oneself that the first two conditions are 
fulfilled. At first sight, however, it appears that there are 
two lengthscales: one associated with the thickness of the 
disordered domains, d (equivalent ly of the interfaces be- 
tween ordered domains), the other with their curvature 
(equivalent to the lengthscale of the ordered domains), 
L. One has to realize that the total interface length / 
is inversely proportional to L. Since the volume occu- 
pied by the interface is conserved we also have d ~ 1//, 
and therefore d ^ L. One thus expects again a growth 
exponent z = 3. 



The above discussion not only applies to Q2 quenches 
(1/2 < Ci < Co-), but also to Q3 quenches {ca- < Ci < Cs). 
A uniform configuration (ccV^o) is then stable, but for 
symmetric initial conditions there are always interfaces 
where = and where it is thus more favorable to have 
c < Ci. In the "ordered" domains c will then increase 
such that we will end up with all the three phases. 



All this is confirmed by our numerical results. In Fig. 6 
we show snapshots of order parameter and concentration 
at late times for q = 0.5 and q = 0.75. In Fig. 7 we 
plot, for the same concentrations, the typical lengthscale 
of ordered domains and disordered domains. They be- 
have similarly. We can also identify two transient eflPects. 
Initially, rapid growth takes place due to the ordering 
process. After that, when the antiphase bounderies are 
wetted, the growth will actually be slower than L ^ t^/^. 
Since the disordered phase has a large surface compared 
to its volume, surface diflFusion will dominate at short 
times, and this is known to lead to a L ^ t^l^ growth 
law [26]. 



4 




C. Ordered regime 



FIG. 6. Snapshots of coarsening in interconnected regime 



(a) Concentration field and (b) order parameter field at 
c = 0.5 at t = 200 (c), (d) id. at c = 0.75 at t = 120 
System size and greyscales as in Fig. 4. 




10 10 10 10 10 10 10 10 10 10 
FIG. 7. Domain growth in interconnected regime (a) Lc vs 

t ,(b) vs t both at c = 0.5, (c) Lc vs t ,(d) vs t both 

at c = 0.75 The local slopes in the insets are calculated as 

explained in Fig. 5. System size 4096 x 4096. 



For Ci > Cs = 1 (Q4 quenches) no spinodal decomposi- 
tion occurs. If initially a domain of the disordered phase 
exists, it disappears. Although the concentration is not 
completely uniform but somewhat lower at the interfaces 
of the ordered domains, as seen in Fig. 9, this does not 
have any influence on the coarsening of the order pa- 
rameter. Thus, the model A exponent z = 2 should be 
observed for q > 1. Again this is confirmed in our sim- 
ulations. We have measured the typical domain size for 
the order parameter around Cg. As can be seen in Fig. 9 
a growth exponent z = 3 is still observed for q slightly 
below Cg, at Q = 0.95. At q = c^, growth is clearly 
faster but the local exponent (inset of Fig. 9c) does not 
reach | in accessible times. The behavior, though, is 
consistent with z = 2. This is most clearly observed 
plotting L againt y/t (Fig. 9b). Increasing q beyond Cg, 
the crossover is now short enough and the z = 2-behavior 
is apparent even in a log-log plot (Fig. 9d). 




c::> 



FIG. 8. Snapshots of coarsening in ordered regime 
(ci — 1.5) at t = 50 (a) concentration field (b) order pa- 
rameter field. System size and greyscale as in Fig. 4. 
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(a) (b) 




10^ 10^ 10^ 10^ 10^ 10^ 

FIG. 9. Growth of typical order parameter domain size 

(a) a = 0.95 (b),(c) a = 1.0 (b): L vs. t^/^ (d) a = 1.5. The 

local slopes in the insets are calculated as explained in Fig. 5. 

The system size is 4096 x 4096. 



IV. DISCUSSION 

A. The microcanonical 0^ model 

The microcanonical 0^ model has lately received re- 
newed attention largely because it offers an interesting 
bridge between statistical mechanics and deterministic 
dynamics [16,15,7]. Its equations of motion derive from 
the well known lattice 0^ Hamiltonian. They can be 
written: 

3 

where the sum is over the 4 nearest neighbors of site i on 
a square lattice. 

As stated before, this model is believed to be in the 
model C universality class since the order parameter 
is coupled to the conserved energy. [1] Indeed, looking at 
snapshots of the order parameter and the local energy, 
defined as: 

one observes that Ei is higher at the interfaces of or- 
dered domains (Fig. 10a) as observed in our model C 
around the coexistence concentration c^. Given the re- 
sults described above, the question of the precise relation- 
ship between the microcanonical 0^ model and model C 
now boils down to whether we are in the interconnected 



((cn + Cs)/2 < Q < Cs) or the ordered regime (q > c^) 
of model C, i.e. whether the interfaces (equivalent ly the 
disordered domains) will thicken or not. 



(a) (b) 




FIG. 10. Snapshots of (a) order parameter and (b) en- 
ergy density E at time t = 15, (c),(d): id. but at t = 75. 



Our observation that, in fact, z = 2 [18] indicates that 
we are in the ordered regime. This implies that the con- 
served field does not undergo spinodal decomposition. 
The snapshots of the order parameter and the local en- 
ergy taken during coarsening (Fig. 10) reveal that the 
width of the "antiphase boundaries" does not grow. In- 
stead, energy flows into the bulk of domains in the form 
of kinetic energy, as observed in [18]. The increase of ki- 
netic energy compensates the energy lost due to the de- 
crease of the total interface length. We believe this con- 
stitutes further support to our result 2: = 2, at least for 
the parameter values studied in [17] and [18]. Whether 
a coexistence regime in the ^^-model exists for other pa- 
rameter values remains, strictly speaking, an open ques- 
tion. In fact, since the system is isolated, the total energy 
and the "temperature" cannot be varied independently. 
Upon increasing the energy, the temperature, which can 
be identified with the kinetic energy, will also increase. 
We therefore only have access to a line in the (c, T)-plane 
of Fig. 1. (Note that increasing the energy corresponds 
to decreasing the concentration.) Investigations of the 
critical behavior of the model, where it is a priori impor- 
tant whether one enters the disordered region from the 
coexistence region or from the ordered region, are under 
way [28]. 
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B. Interfacial properties 



A few years ago, Somoza and Sagui have investigated 
the effect of interfacial properties on the morphology of 
domains in a model C system [22]. They came to the 
conclusion that apart from a "complete wetting" regime 
where two "ordered" domains {ip = ±1) are always sep- 
arated by a disordered layer {ip = 0), there exist also 
a "partial wetting" and even a "partial drying" regime 
where ordered domains of opposite sign can be in direct 
contact. These regimes were observed by varying the dif- 
fusion constants Kc and keeping Ax and At fixed. 
Thus diffusion of one quantity was enhanced with respect 
to the other. 

We first repeated the simulations of Somoza and Sagui 
using Ax = 1. (Fig. 11a). The ordered domains are in 
direct contact but the corresponding interfaces are not 
smooth, an indication that the equations are not well re- 
solved. Next, we performed the same quench but at a 
higher numerical resolution {Ax = 0.3). As shown in 
Fig. lib, no direct contact between ordered domains can 
be observed. The argument put forward in [22] is that in 
the partial drying regime the surface tension is lower be- 
tween the two ordered phases than between the ordered 
and disordered phases. This surface tension was mea- 
sured by calculating the free energy of a numerically ob- 
tained interface profile [27]. It must be noted, however, 
that this procedure is only correct when the obtained 
profile is the equilibrium profile. An order-order inter- 
face for instance will desintegrate into two disorder-order 
interfaces when the surface tension of the latter is lower. 
The speed at which this happens depends on the diffusion 
constants, since some amount of the conserved quantity 
has to diffuse away from the interface. One could thus 
have the impression that an order-order interface is sta- 
ble, when in fact it is very slowly disintegrating. 



(a) (b) 




(c) (d) 




FIG. 11. Snapshots of partial wetting regime = 0.6, 
= 1/Kc with Ax = 1.0, At = 0.01. (a) Concentration 
field, (b) Order parameter (c),(d), Same system but better 
resolved: A = 0.3, At = 0.005. System size and greyscales 
are in same as in Fig. 4. 



V. CONCLUSION 



We have provided a comprehensive account of do- 
main growth in two-dimensional model C systems, uni- 
fying and sometimes correcting partial results present 
in the literature. Quenches into the coexistence region 
{cn < Ci < Cs) all lead to 2; = 3 domain growth governed 
by the conserved field, independent of the morphology of 
domains, but subjected to the strength of initial fluctua- 
tions when Ci < ci . Quenches into the ordered region are 
asymptotically dominated by the ordering process, lead- 
ing to z = 2 growth after a crossover from slower z = 3 
coarsening. 

We also clarified the status of the microcanonical 0^ 
model and showed that in this case sub-critical domain 
growth correspond to quenches into the ordered region of 
model C. 

The picture that thus emerges is that domain growth 
in model C systems below criticality is governed by the 
growth exponent of either model A or model B. At crit- 
icality, on the other hand, one may expect exponents 
different from model A or B. The natural question that 
then arises is how the exponents depend on on the initial 
concentration. 

Another issue open for further investigation is the be- 
havior of the so-called persistence probability p, defined 
as the probability for a spin to have remained in its initial 
phase from to up to time t (see [29] for a review). In the 
case of model A and B this quantity decays algebraically 
with an exponent 0, estimated to be Oa ~ 0.21 [30] and 
Ob ~ 0.25 [31]. The universality of these exponents is an 
ongoing debate. Measurements of the persistence expo- 
nent in model C, especially in the interconnected regime 
where one would expect the model B value, might help 
to settle this issue. 
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